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EXECUTIVE SUMMARY 

As pirt of the Ottawa River nuclear spin contingency planning studlfis, a model 
which will aid in predicting the arrival times, durations of passagi, and 
concintratlon distributions of the spnied radionuclide at the water supply 
intakes and other strategic locations has been developed. The model utilizes an 
analytical solution for the two-dimensional finite time case, the derivation of 
which Is based on the stream tube concept wherein the lateral co-ordinate Is 
represinted by the cymulatlve discharge. An approximate Integration procedure, 
known as the Laplace method, has been utilized in the derivations. The 
analytical solution applicable to pipe outfans inclydes the effects of 
side-waTl ref lection of the material from channel banks. Expressions have been 
presented for the expected arrival times and durations of passage of spilled 
material cloud at strategic locations. 

A model with the acronym, TWODIFIN, has been developed to perform the 
computations. The miodel accounts for longitudinal variabilities of channel 
geometric and hydraulic parameters (viz, width, depth, velocity, longitudinal 
dispersion coefficient, and lateral diffusion factor) and radlonycllde decay 
rates. This Is achieved by dividing the channel study segment Into a number of 
reaches such that the transport process parameters are fairly constant within 
each reach and by ut ill zing the moving average values of the parameters in the 
computations. The concentrations and flow rates of the spill discharge are 
specified as step Inputs occurring over a known number of time intervals. The 
other Inputs to the program include the channel hydraulic parameters, 
longitudinal dispirsion coefficients, lateral diffusion factors, and 
ridionuclide decay rates for successive reaches of river channel below the 
outfall. The output of the model includes lateral concentration distributions 
of radionuclide material as a function of time at specified locations below the 
outfall , The times of prediction are generated from the calculated arrival 
times and durations of passage. 

An example based on the data gathered in the Ottawa River following a spill of 
tritium from the Rolphton nuclear power demonstration (NOP) facility during 1981 
illustrates the input parameters and output summary tables of the TWODIFIH 
model. It should be noted that this example is an illustration of the 
computational aspects; it is not intended for a validation of the model. 



Tht importance of the convective period 1n the Ottawa River below the Rolphton 
NPD dlschtrfe has been assessed with the aid of the tritium 
concentratlon-versus-tiiiie distributions monitored at Deep River, Petawawa and 
Pembroke during August - September 1981, following the trltiui spill from the 
Rolphton NPD in August 1981. The convective period Is found to dominate for 
about 16-20 km below Rolphton as indicated by a skewed concentration 
distribution at Deep River; however, the digrte of skewness relative to that 
for a normal distribution (skewness ^ 0) Is found to be small. Thus, it is 
concluded that the effect of convective ptrlod can be neglected In the modelling 
studies. It is recommended that this aspect be further examined with the aid of 
data from the recent field studies (condycted during September 1984) and 
modening studies. 

The relation bttwein the convective period and lateral mixing zones has been 
examined. A detailed assissment could not be made because of a lack of relevant 
Information (either theoretical or field study data); it 1s recommended to 
further evaluate this aspect with the aid of 2-0 model and the field data 
gathered In the recent field studies. 

Methods for Unking near-field and far-field models have been presented. A 
schematic representation shows the various zones in the near-field and far-field 
regions, and the models applicable to these zones. The methods for linking of 
models presented herein will permit a continuity in applying the models for the 
near-field and far-field regions. 



fl 



L*'ilabo^tion du plan d^xttgenoe si 'Cas. de. deveisaneiit riiclmire ^dans la 
riviere des Outacuais nam a amenes a nettre au point mi modele ^li aidera a 

prMire l*.arrivte des. radionucleid^ aux prises d'^eau et a d'auto^ points. 
strategiqu^, leur 'durte de transit et la disteibutlai, f3m Imam 
'^ncentrations. 1b model© utilise une solution .analytiqu.e poor desire 
1' incident bidimeisionnel a teips fini, iSolirtiOTi ri^xMant mjx le oono^ft, d^. 
ti&es d'eosulQiiadt dans lequel la coorfonnie latirale ^^it n^riseTtie pa^ le 
devea^aneiit oiamlatif . Pour parvaiir a»c dfei'^ies, nous nous semn^ servis 
d'une ■t«dmiiqi2e d'integiation a^roKinatlve oamue sous, le ncm de ■iiitliode de 
laplace. La :BQlution .analytique qui s'.^pli^ie: mm points de dev^ssKient 
tiait corpte d^ nucleides. ^li sont, reflets d^. berg^. du dienal. D^. 
esipress.iQns ont ete foirnulees ■pour' I'a^ivte 'estiinative Ai :nuage ai d^. 
ai^r>its cles. a,ins.i igue sa, 'duree. probable. 

Un nodele portant l"a.eronp[ie 'TOODIF^I a ite el^abori pour effectaaM' 1^. 
^iloils. 'Ce. inodele rend ooipte d^ variability. longitiMiiial^ des. paranitres. 
georaitri^^ et hydrauliqu^ du Aenal Cc'''^5t-A"dire la largeur, la 
profondeur, la vitasse, le ooefficiait de diversion lo^itiidijBle et le 
facteuT' de dif:fus.ion laterale) ainsi que d^ taux de dfeintigration d^. 
.radionucleides... II y .arrive en. divisant la partie ^du Aenal k I'etajde ei, un, 
certain nanbre de tron^ns dans l^^els lee pazOTetrM dH d^laoement soTt 
.assez oonstants,, puis ai utilisant les. valairB ma^&wmm du .d^laoemsnt dans 
les calojls. Lbs ooToaitratioi^ et les taux d'eaoulaiient des si^bstano^ 
dev^Bies. ;sont oonsideres oanroe des. eitm^ qpi se produisent a un, :noiibre 
cGimu d''interwal„l^.. Ibb .autr^. 'Sntre^ sont 1^ parmnitr^. hydrauliqu^ du 
ctieaial,, les. coefficients de 'dispersion longitoainale, 1^ facteurs •de 
diffusion lat&ale et 1^ taiK de d^infcegratlon des radionucleid^ pour 1^ 
toongoiK suooessifs du, '(Aenal, 'de la rivifere ^en mmL &ii :point. de ^EvaraemaTt, 
0iant aux sortie, ell^ oajpjDaineift la distributiOTi d<^ oonafflTtrations de 
radionucliides, ei f auction du tenps, a d^ witooits precis mn. aval du point 
de deverBement. L^ pnidictiais ta^oBBll^ provieraiaTt du moraHTt d'a^ivie 
et de ,1a durte: de 'trans.,it, dija '^OLoilto.,. 



'Uh exQiple 'qui s'lJi^ire d^ domn^^ ZBCueilli^ dans la rivitoe d^ 
Outooiais pff suite du dw^seniait de tritiiii suryerm ai 1981 a la oentrale 
nuclea,ire P'ilote de lolphton illustre bioi les ^raflmetr^ des 'fflitiiies at 1^ 
tableaux ricapi-bilatifs d^ sortie Ai lodele TOODIFTO. II faut faire 
ETemartper que cet esfla^la a pour seul but d'illustper les calculs, et rion de 
verifier' le nodele, 

L' importanoe de la periode de cOTwaction dans la riviere des C3iitacuais en aval 
du deveirseiiaTt, ;syrvenu i la oaitrale de Bolpliton a ete waluie au wa^mi des 
■OQTipentrationB de ■tritium,, mi. fonction du, taips, qui ont ete, mesuri^ a Dei^' 
Riv^-, a Betairawa et a Penbroke en aout et &rk s^taitore 1981, apris 
I'incideit, qui s'est proAiit &n aoat de la niiie annte* la periode de 
convection, avone-nous constate, doraine sur mm distance ^^iroaciiiBtive de 16 a 
20 ki,l,CTOtr^ en, aval, da Bolphton, oanone le montre, la 'distr'iJ3Ut,ion a^metriqua 
des ODncentrations a Deep :Rivari t«mitefois,- la degre d'a^mitrie ^est idnime si 
on le ccMpare, a celui d'''une dietribution nomale (val,eijr 0') ,. Ainsi, nous 
'Concluons qu'il n'^it ^s naa^^aira de toiir cooopte des^ effets^ de la, p^iode 
dfe oonvection dans la reaHsation des iiodeles. A cat egazd, nous r^^cramandons 
■^e Get a^ect so,it acmdni plus avant ai sa s^ryant 'd^ donnees. provenant d^^ 
etudes rialisees :sur le tenrain en s^taribre 1984 at des etudes de. 
modelisation. 

:La relation ^^li ^ciste, ®itre: la periode de ocjfTvection et les. zoo^ de melange 
,laterales. a eti: etiMife:, Une evaluation detaillee s'^est averte jiipoes.ible a 
^use du manque de rKTseigr»iiaTits, soit theoriques, soit pratiques; nous 
pr^osons done ^e ^^t elanait soit approfoaidi au n^ai du modale 
bidimeTsionnel et des donnas r^^^eillies au cours des etud^ sur le teixam. 

Ijes. moyens. de lier 1^ models 'du Aanp prodhe 'et 1^ models du 'drianp 'iloipie 
ont 6te prteentte, la r^rteentation sdiemati<pe Indique 1^ div^^^ zon^ 
qj± sa troyvoit dans 1^ deuac ^lemps, ainsi qua 1^ nodeles qui s'y 
q^li^iint, o^ nifto^s de liai^n d^ modal^ periiettoait d •assure' la 

oraTtinuite d' application de oemc-ci daj^ le ^is at d^ dia^is pi?odi^ et ^ss 
Amps iloi^is. 
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I INTRODUCTION 

1,1 BacikgirO'Uind 

Contingency planning studies undertaken by the Ontario Ministry of the Environment 
(NOE) for the management of accidental spill s or controlled release of 
radionuclides into the Ottawa River from nuclear power pTants include divelopment 
of mathematical models having the abnity to predict the concentrations and 
durations of passage of such radionuclide releases at municipal water intakes and 
other strategic locations. As a part of this study, methods for estimating the 
longitydlnal dispersion coefficients as well as an analytical equation for the 
one-dimensional finite-time time release case, 1.e. transient 1-D model , have been 
presented in a previous report (Gore & Storrle, Ltd., 1983). Some of the aspects 
requiring further studies which were identified in that report Inciude developiient 
of a far-fiild f1n1te-tiie 2-D model, methods for Unking near and far-field 
models, assessing the importance of the convective period and evaluation of the 
relation between the convective period zone and lateral mixing zone. 

A finite-time 2-D model is required because of the fact that the discharge at the 
source is either from a pipe or a diffuser outfall which usually results In lateral 
concentration gradients being set up for a considerable downstream distance. For 
example, the distance required for complete lateral mixing downstream from the 
Rolphton NPD is esttroated to be nearly 43 km; on the other hand, vertical mlxinf 
1s esttmated to take place in less than 1 km (Gore I Storrie, Ltd., 1983). It was 
also estimated that the 1-D bulk diffusion equations, valid in the diffusive 
period, are applicable after a convective period which dominates for about 50 km 
below the Rolphton NPD, The length of the convective period zone was estimated 
from an equation developed by Fischer (1967); however, this was not assessed with 
the aid of concentration distributions. These estimiates also indicate that the 
zones of convective period and lateral mixing extend for aliiost the same length of 
the river, which points out a need to assess their relative significance In the 
context of modelling studies. Nethods for linking near-field and far-field models 
are also needed to establish continuity among near-field, far-field 2-0 and 
far-field 1-D zones in modelling studies. 



1.2 Obiectives 

The objictlves of this Investigation are as follows j 

1. To develop a far-fielcl 2-D finite-time model, which Includes longitudinal 
advection, dispersion and decay, and lateral dispersion. 

2. To present liethods for Unking near-field and far-field lodels. 

3. To assess the importance of the "convectlve period" in the Ottawa River 
below the Rolphton NPD. 

4. To evaluate the relationship between the convective piriod zone and the 
lateral mixing zone. 



U FAR-FIELD 2-D TRANSIENT MODEL 

f.1 Theoretical Aspects 
2.1.1 Differential Equation 

The basic differential equation describing the conservation of mass for the 2-D 
transient case can be obtained by following the procedures presented by Yotsukura 
and Sayre (1976). The differential equation, based on the "streain tube" concept 
wherein the transverse co-ordinate is represented by the cumulative discharge (q) 
can be written In the following fori: 

3c u, 9c e^ 3^c yD„ d^c 
9t m 3x m^ 3x^ m^ 3q^ 

XX A 

where x = longitudinal co-ordinate (m) 

q s transverse co-ordinate (m /s) 

t - time (s) 

Ux - velocity in the x-direction (m/s) 
nix ' metric coefficient (or scaling factor) associated with 

X-axis 
ex - longitudinal dispersion coefficient (m/s) 
Dy = transverse diffusion factor (m^/s ) 
K^ - first-order decay rate coefficient (1/s) 

c = concentration at a point (x, q) at time (t) 

Equation 1 is subject to the following assumptions: 

1. The density of effluent and receiving water are the same. 

2. The concentration distributions in the far-field region are not affected by 
the near-field processes such as the jet-effect. 

3. The vertical concentration distributions have attained uniformity just below 
the outfall . 

4. The decay of pollutant follows the first-order kinetic law. 



The independent variable (q) and the diffusion factor (Dy) appearing In 
Equation 1 are defined by: 



q = J my Uy hy*dy (2) 



Dy = mx By Uy h^ " ff ^^ ^^ "^ "^y"^^ -^- 



whereln the overbar indicates the average of the product terii 
(m^ ey Uy hy); my is a metric coefficient associated with the 
y-ax1s; hy and Uy denote ^ in order, the local depth and local 
depth-averaged velocity 1n the x-d1rect1on at a lateral distance (y) w.r.t a 
refirence bank of channel; ey = transverse dispersion coefficient; and 
Q = discharge 1n the river channel Just below the source. The metric coefficients, 
fflx and ffly, account for the convergence or divergence of stream surfaces 
along the respective axes, x and y. The diffusion factor (Dy) can also be 
expressed by (Gowda, 1980; 19844; 1984b) ^ 



Dy - mx f ey y ti^ - Q-6/b (4) 

where h and u denote, in order, the cross-sectional mean values of depth and 
velocity; b = channel width; f = shape-velocity factor; 
B = nondlmensional diffusion factor; and ix - average value of mx- 

Net hods for estimating ex» Dy, B, f » Wx and ™y l^^^e been 

described in various publications (Yotsukura and Cobb, 1972; Beltaos, 1979; 

Gowda, 1980; Gore I Storrle Ltd., 1983). It should be noted that the expressions 

for B and Dy incorporate f , mx and my. A summary of the values 

of B for natural streams can be found In other publications (Gowda, 1981; 

1984a). The decay rate coefficient (K^j) at a known temperature (Ti°C) 

for a given pollutant may be obtained from field study data or a literature review. 

The rate coefficient at a design temperature (T2®C) is calculated from the 

modified van't Hoff-Arrhenlus relationship: 



T2 ~ Ti 

Kt2 = Kjl i " fi) 

in w^hlch Kxi ind %J2 denote the ratt cotfflcients at T|° and 
Tg" C, respectively; and 8 is the temperature correction factor. 

2.1.2 Analvticsl Solution: 

For the sakt of convenience in obtaining analytical solutions, Equation 1 is 
rewritten as follows: 

Sc 9c 3^c 3^c 

_ + u — = E — + E - K .C (6) 

it ax" ^ dK^ ^ dq^ - 

where 

u e uD„ 

Although mx occurs 1n these equations, it can be assuned to have a unit value. 
This approximation 1s not likely to have a major effect on the computations. 

An analytical solytlon of Equation i for the case of an instantaneous pollutant 
source located at x « and q = qg (with respect to origin at a reference bank 
of channel) Is given by: 



Cj = -p^^ exp 



(x - ut)^ (q. - q)^ 

• = - =^= - K.t } (8) 

4E^;t 4Eqt 



1 



In which q^ ^ cumulative discharge between the reference bank and the source; 
Wi = total mass of material discharged instantaneously; and 

Ci = concentration at a point (x, q) at time (t) due to Uf, 

The terms occyrring inside the square brackets in Equation 8 can be simp Tif led as 
outlined by Glover (1964) and Wnek and Fochtman (1972) to obtain a relationship: 



^1 = 



_ W'|U.exp (uix/2e^) exp ( - | - nt) 



(9) 



2wJO 



X q 



.hare C J jl , i%:j}^) ^ and n = j^ , K, j (^q. 



,4E 4E I (4E 



Thi concentration distribution due to the release of a neutral density material at 
a rate (w) at the point x - and q - during a finite time (dt') extending from 
t' « tg to t' = tn can be obtained by Integrating Equation 9 with respect 
to time. Thus, for the finlte-tiiie case, we have: 



exp (ux/2ej f " w.df ( / ^ \ 1 



(11) 



The rate at which the material is released (w) may vary over the duration of 

injection (tn - to). In such cases, the duration (tp = tg) can be 

divided Into a number (N) of infinitesimally small time intervals such that within 

each infinitesiiial interval (tj - ti), the injection rate (Aw) Is 

fairly constant. Then, the concentration (Ac) at a point (x, q) at a given 

t1«e (t) due to the constant rate injection during (tj - t^) is given by: 



,„ u.exp (ux/2e ) fH df f / ? \ -1 

^'' " ^^- (Ay) I =— ^. . exp . . Ti^ (t ^ V] (12: 

t. 



Thi total concentration (c) at the point (x, q) at the time (t) due to the iiterial 
released over the entire duration (tn-to) - N (tj - t^) can be 
obained by summing the Individual Ac values. Thus: 



N 



x^O 






The time integral appearing in Equation 13 can be evaluated approximately by a 
procedure known as the Laplace lethod (Erdelyi, 1956), the details of which can be 
found in Appendix A. This method has been utniied to solve 2-0 and 3-D equations 
by Henry and Foree (1979) and Sagar (1982). However, the limits of integration in 
these two studies have been taken to be ti ^ and tj = t (I.e., from to 
t), the upper limit of which coincides with the time (t) at which the concentration 
at a point (x, q) is calculated. Wnek and Fochtman (1972) have presented 
analytical solutions for the 3-D case without resorting to the approximate Laplace 
method, but by a direct integration procedyre using the integration limits, to t. 
Obviously, the use of Integration limits to t is inconsistent with the solution 
being sought. Concentrations calcylated from those analytical solutions are either 
=ve or unrealistlcally large* Medina and Buzun (1981) obtained analytical 
solutions for the 1-0 case by direct integration using the Hmrts to ti. 
Such a procedure was also the basis of development of analytical solution for the 
finite-time 1=D case in the Ottawa River Study (Gore I Storrle Ltd., 1983) i 
however, it was pointed out to evaluate the need for changing the integration 
limits to ti to t2. These solutions (i.e., with integration limits to 
ti) are also found to predict unreasonable concentration values. Based on 
these analyses, the integration limits tj to tj are utiTized 1n the 
formulations presented in this report. 

Evaluation of the integral in Equation 13 by the Laplace method (Section A.l, 
Appendix A) will result in the following expression for the concentration 

distribution: 




^^^^'l-dF-/^) --(JT-/^)11 



114a.} 



where T, = (t - t.) ind Tj =• |t - tj 



iind T, =• (t - tA. Iiail 



2.1.3 Analytical Solution with Side-Wall Reflection Effects: 



Equation 14 does not include the effects of side-wan reflection of the material 
from the channel banks. Following Yotsukura and Cobb (1972), let us consider that 
the material Is released as a point (or vertical line) source at a location x = 
and q - q^. The solution which accounts for the side-wall reflection of the 
material is then obtained from the method of images. For the sake of convenience 
in writing the resulting solution in a compact form, the following terms are 
defined (assuming m^ - 1): 



(^ ■ -.) ■ 



iQ^u 



P - q/Q ; 



%/Q 



fi - fa(2n + Pg -p)'7'4i0 u 

fi - b(2n + Pg +p)|2/40 y 

f3 = b(2ni - p,g -p)V46 u 

'fu = b(2in - p,^, +p)V4B u 



'^ b 

(15b) 

(15c) 
(15d) 

(15e) 
(15f) 



(15a) 



where n denotes the number of images. Then, the analytical solution which includes 
the side-wall reflection effects can be written as follows: 



b u 



4Q ^ TTiBe^ 



N 



:n' 






+ E 



n=l 



4^)- ' Fi 



a + f 



{^^a * t)}'^'^ 



\-(J^f-IT.)-"(J^^-l^)\ 



"p(5^)-2 J.(a.,jj . p-^ . / n^ r—\\ 
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In solving Equation 16, Ji or Tj could be zero or negitive 1n some 
computational stips depending on the values of t, t^ and tj under 
consideration. When Tj (or Tj) 1s zero, the argument of the error function 
becomes infinity; and hence, the vtlut of the corresponding error function will be 
[iff(ffl)] = 1.0. However, when Ti and/or Tj attain negative values, the 
arguments of the error functions become imaginary; for such computational stepSs 
the concentration values cannot be evaluated. This is consistent with tht fact 
that the inputs at the times t^ < t and tj < t do not contribute to the 
concentration distribution, c (x, q, t). 

There coy Id be an upstream movement of discharged material when the diffusive 
transport is greater than the longitudinal advective transport. For such 
conditions. Equations 8, 14 and 16 are applicable to negative as well as positive 
values of the longitudinal distance (x) with respect to the origin at the source. 
The distance affected by the upstream movement (i.e., negative values of x) Is 
dependent on the values of e^ and u, 

2.1.4 Times of Passage; 

A knowledge of the expected time of arrival and duration of passage of cloud at a 
given longitudinal distance (x) from the source 1s necessary to predict 
concentrations using Equation 16. These can be estimated through the use of the 
relationships applicable to an instantaneous source 1n conjunction with the known 
duration of the finite-time source. The passage times of cloud in the 2-D 
instantaneous release case can be determined from the following relationship (see 
Section A. 2 in Appendix A): 



t . -iyj ux + le^ In Cl/f"} t I Jt ,u X e^ In (l/r)]i+ [ e^ In fl/T') ], "^ - | e^ ub (q^ - qVim'l \ 



(17)1 



where f ' is a specified relative concentration. It denotes the approximate 
ratio of the desired minlmym concentration to the peak concentration at a given 
location. The selection of Its value is dependent on the desired accuracy , the 
general range being 0,1 - 0.01. 



Iquatlon 17 gives the passage times as a function of the cymulitlvi discharge (q) 
at desired transects below the source. For computational convenience, the times 
obtained with qj = c| (i.e., the maxiwum passagi tiraes) at a given location, 
will be considered for all values of q. It should be noted that Equation 17 
obtained for the 2-D case will be the same as that for the 1-D case when q ^ q^ 
(see Heathcote, 1982). 

Froi Equation 17, two values (ta and %) are obtained with the -ve and +ve 
signs, respectively. The value for t^ corresponds to the arrival time; indl 
the difference {% - tg) gives the duration of passage of the cloud for the 
Instantaneous case. The expected arrival time of cloud for the finite-time case 
will also be equal to tj. However, the duration of passagi (tp) should be 
modified to account for the duration of the finite-time input; thus, the resulting 
expression is 

tp = tb - ta + tf (18) 

In which tf « duration of input from the finite-time source. The procedures 
utilized in computer program development to determine the number of time steps and 
the times for predicting the concentrations at a given transect are described in 
detail In Section 2,2, 

2.1.5 Reach-Dependenicv of Parameters: 

In natural streams and rivers, the channel widths, depths and velocities, 
longitudinal dispersion coefficients, lateral diffusion factors and pollutant decay 
rates tend to vary 1n the longitudinal direction. In such cases, moving average 
values of the parameters are utnized in the coiputations (Gowda, 1980; 1984b), 
The pirtinent relationships are presented below. 

The moving average values of width, i-j, depth, H-f, and velocity, Ui, 
applicable to the channel length, x^, are determined from 

1 4 B^x^ 

li 



The relationships for moving average values of longitudinal dispersion coefficient 
(6x1 ) *"<' transverse diffusion factor (Si) are 



i 
= L ^ h'vi (^ - xui) (20) 



^xi = r. jii ' -^ ^ J " J-1 



'r- 7^ A le'j (Xj-Xj.i)l (21) 



where e'xj and 8'j denote the values for the jth reach of length 

The moving average values of the decay rate coefficients (Ki) are obtained froii 

J<1 - 2 (K'j Ttj)/Tti (22) 

wheri K'j Is the decay rate in the jth reach; Ttj is the travel time In 
the jth reach of length (xj - Xj-i) and T^i 1s the total travel time 
in the length x|, extendi ng from the source to the transect, i* 

.2.1.6 Estimiation ol the Nurnber of Ilm:aiges: 

An examination of Equation 16 Indicates that the virious summation terms Involving 
the images can be teriiinatid by setting: 



expl/— \- 2 ^T](a + f,^) I = 6 ; m = 1 to 4 (23) 



where 6 1s a positive number to be selected such that the value of each nth 
term is small enough to cause a negligible effect on c(x, q, t). This procedure is 
similar to the one utilized in mixing zone models (Gowda, 1980). Taking logarithms 
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on both sidis to base e and simplifying, one obtains quadratic equations of n. The 
number of Imagss for the four tens involving f|p (m = 1 to 4) are then 
given by the fonowing solutions of the qytdratlc equations: 



f^: n ^ 0.5 (p - ps) ± X 
f2: n = -0.5 (p + Ps) ± x 



fj- n - 0.5 (p + Ps) ± X 



¥4: n - -0,5 (p - Ps) ± X 



(24a) 
(24b) 
(24c) 
(24d) 



where A 



m 



I i 



[ (ux/2e ) - In (1/6) ] 2 



^h 



- 4a 



(24e) 



subject to the condltiom that [ (ux/2e, ) - In (1/6) ] ^ ; 4rp (24f) 



1.e, 6 > exp 



m- ^ 



(24g) 



Generally, the valyes of 6 vary approxiiiatily betwtin 1.0 and 10"^ 

for sman and large values of x, respectively. 

It is recomiiended to determine the value of S for a given set of conditions as 

outlined tlsewhere (Gowda, 1980), 

FroB Equations 24 a - d, two values of n are obtained for each case: A lower bound 
n| with the negative sign and a higher bound 02 with the positive sign. 
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The total number of images, nt, for each case is given by the sum of the 
absolute values of nj and n2. The sunmatlons, n ~ to ® and n = 1 
to ®, appearing in Equation 16, are then caried out froii n * to n^ and 
n = 1 to nit, respectively. Generally, the maximum number of images requirid is 
nine (Gowda, 1980). 

2.2 Computer Program Development 
2.2.1 General Description 

A computer program for the TW0-DImens1onal FINite-time model with the acronym 
TWODIFIN has been developed to perform the computations according to Equations 
16-24. A listing of the program, written in FORTRAN language, 1s presented in 
Appendix B. A list of the input/oytput variables and an example showing the input 
and output of the program are also presented in the appendix. The salient features 
of the computer program TWODIFIN are presented below. 

The inputs to the program inclyde streamflow rate and the corresponding width, 
depth and velocity, as well as nondimensional transverse diffusion factor and 
pollutant decay rate for each reach of the river segment under consideration, and 
the longitudinal distances from the source (i.e., oytfall) to each transect. The 
design terapirature and streamflow values are also specified. The widths, depths 
and velocities for each reach are scaled up to the design streamflow rate using the 
Loepold-Maddock equations, while the decay rates are adjusted to the design 
temperature according to the modified van't Hoff-Arrhenlus relationship (Equation 
5). Then, the moving average values of these parameters are computed. 

The longitydinal dispersion coefficients (ix) for each reach can be directly 
Input to the program by specifying the flag parameter, FLAGEX = 1. Alternatively, 
if FLAGEX = 2, the ©x values are computed from (Bansal, 1971; Gore & Storrle 
Ltd., 1983) 

fij^ - [(7.05)(10^) u R / R^'"^^ 1 mVs (25) 

where R^ - ( u fi/v) is the Reynolds number, v being the kinematic 
viscosity of water (m /%} at the known temperature (Ti° C), 
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The viscosity at a design temperature (T2° C) is calculated from 
(Rich, 1973) 



The discharges, pollutant concentrations and associated times at which the inputs 
enter the river channel, as well as the location of pipe outfall in the river 
channel (as a fraction of the river discharge), are also inclyded in the input 
data. Other data inputs to the prograii are stated 1n Appendix B. 

The duration of passage (tp), estimated from Equations 17 and 18 as outlined 
earlier, Is utilized to determine the nuiiber of time-steps (nt) as follows: 

If tp i 10, nt - 20 

If 10 < tp i 40, nt = 40 

If tp > 40, n| = 50 



The numbers of time steps for the rising and falling portions of the 

concentration-time curves (nn and nt2t respectively), are calculated from 

the empirical relationships nti = (0.4 n^) and nt2 - (0.6 n^). 

This empirical procedure is designed to permit concentration predictions at 

relatively closer intervals on the rising portion of the curve compared with the 

falling portion which tends to display an elongated tail. 

The approximate time (tf) to the concentration peak at i transect distant (x) 
from the outfan 1s calculated from 

tr = [(0.8 x/U) - ta] hours (27) 

where U = velocity (m/hr) and tj = arrival time (hour). 

Then, the time-Intervals for the rising and falling portions of the curve 
{At| and 4t2, respectively), are obtained from 

Ati = tr/(nti = l)t (28a) 

m 



and At2 - (tp - tr)/(nt2 - 1) (28b) 

The ttroes for predicting the concentrations are then determined from the 
expressions 

t = [ta +■ m (k - 1)]. for 1 I k I ntl, (29) 

t « [t|. + At2 (k - ntl)], for ntl < •< i nt2. (30) 



The output of the TWODIFrN prograni includes the distance from the source to the 
transect under consideration, the associated moving average values of the process 
parameters and the lateral concentration distributions at various passage times at 

the transect. 

It should be noted that the algorithms of TWODIFIN are valid for positive values of 
X, although the various analytical solutions (i.e,, Equations 8, 14 and 16) are 
applicable to negative values of x as well. 

2.2.2 Exannple 

In order to niustrate the data inputs for the TWODIFIN model and the resulting 
output , an example based on the data gathered following the tritium spill from the 
RoTphton NPD during August 17 - 25, 1981, is considered. The input data for this 
example and the computational resylts are presented 1n Appendix B.. The salient 
aspects of this Illustrative example are discussed below. 

The number of transects considered is five (Inclydrng x « 0). These Include the 
locations at Rolphton, Deep River, Chalk River, Petawawt and Pembroke (see Figure 1 
1n Chapter III). The channel hydrauTic data are known at a flowrate of 298.7 
mVs. For the tritium spill conditions of August 1981, the streamflow rate is 
493.1 mVs. The exponents for scale-yp of widths, depths and velocities are 
taken equal to 0.375, 0.375 and 0.25, respectively, for each reach. The 
longitudinal dispersion coefficients are computed Internally in the program, using 
Equation 25. The value of nondimenslonal diffusion factor for each reach is 
assumed to be 0.001. The decay rate of tritium is taken to be zero. The relative 
concentration (¥') is set at 0.1. The flow rate and associated tritium 
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concentfatlori values at various tines wers obtained from the concentration versus 
time plot monitored at the Rolphton NPD (sit Figure 2 in Chapter III). 

The sequence of Input data for the prograni 1s described 1n Section B.3, Appendix B, 
The Input data for this niustratlve exaniple are given in Section B.4 of the 
appendix. 

The locations at which the tritium concentrations are computed by the TyODIFIN 
■odel Include Deep River, Chalk Rrver, Petawawa and Pembroke. The output tables 
for these four locations are presented In Section B.4, Appendix B. These tables 
show lateral concentration distributions at q ^ to 1 (with an interval of 0.2) 
for various times. The moving average values of the process parameters are also 
rnclyded In the tables. 
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Ill AIMALYSES RELATED TO CONVECTIVE PERIOD 

AND LINKING OF NEAR-FIELD AND FAR-FIELD MODELS 

E'l Assessment Qf the Importance of Convective Period 

The theoretical distance downstreaiB of the Ralphton HPD discharge within which 
convective period 1s likely to be dominant, has been estlmatid to be 53.75 km (Gore 
it Storrle, Ltd., 1983). This estimate was obtained from a relationship developed 
by Fischer (1967) based on laboratory flume studies. The convective period is 
characterized by highly skewed concintrat1on=v§rsus=t1me distributions. NcQulvey 
and Keefer (1976) have developed a model termed ""convective model of longltydinal 
dispersion" which can be utilized to predict concentration distributions in the 
convective period zone. IHIowever, the need for applying the convective model should 
be determined with the aid of measured concentratlon-versus-timie distributlo'ns In 
the study segment of the river channel. This is necessary since the length within 
which the convective period 1s dominant in natural rivers could differ from that 
predicted by Fischer's equation depending on the presence of meanders, riffles and 
other conditions that enhance the mixing processes. 

In order to assess the importance of the convective period in the Ottawa River 
below the Rolphton NPD discharge, data collected as a result of an accidental 
release of trit1um-contam1nated water from the Rolphton NPD plant during August 
1981, will be utilized. The discharge Involved a total of 3786 Ci (1 Ci « 1 curie 
- 27 X 10^° disintegrations per second) over a duration of about nine days 
(August 17-25, 1981). During the ensuing period, tritium was fflonitored at the 
municipal water supply intakes located at Deep River, Petawawa and Pembroke, as 
well as at Ottawa. Figure 1 shows the study arei including the Rolphton NPD plant 
and the monitoring locations. 

The amount of tritium discharged from the Rolphton NPD plant varied with time as 
shown In Figure 2. Plots of tritium concentration versus time, measured at the 
three water intake locations and at Ottawa, have been presented In Figuri 3. 
(Note: In Figure 3, the tritium concentrations are in BecquereTs/L where 1 
Becquerel/L = 27 x 10"^^ Ci/l). This figure, prepared by the Radiation 
Protection Bureau, Health and Welfare Canada (1983), shows the 
concentration-versus-tlme plots as step-functions. In order to aid in the 
identification of any skewnesses In the distributions, smooth curves (shown by 
broken lines) have been drawn on the plots in Figure 3. An examination of these 
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sioothed curves indicates the presence of some degree of skewness in the trltluin 
dlstrlbutlQn at the first monitoring site at Dtip River (located about 16 km below 
the Rolphton NPD). The distributions at Petawawa and Pimbroke {located about 40 
and 60 kin respectively, below the Rolphton NPD) do not exhibit a clearly 
discernible skewness. Therefore, 1t is reasonable to conclude that the length of 
the convective period zone In the Ottawa River below the Rolphton NPD discharge 
would be in the 16-20 km range. 

A further examination of the detailed tritium concentration-versus-tlrae plot at the 
Deep River monitoring site, presented in Figure 4, indicates the foil owing 

characteristics: 

coefficient of skew = 0.18 
kurtosis = 2,01 

A comparison of these two parameters with the corresponding values for a Gaussian 
distribution (skewness = 0.0; kurtosis = 3) indicates that the observed tritium 
distribution at Deep River can be approximated by a Gaussian distribution for 
modening purposes. This will permit the utilization of convection-diffusion 
models without the need for a convective model. However, it 1s recommended to 
further evaluate the implications of this approximation through a comparison of 
observed and predicted durations and time of passage, and concentration 
distributions for the Deep River site, using the data collected in the study aria. 

3.2 Relation Between Convective Period and Lateral Mixing Zones 

The longitudinal distance below the Rolphton NPD discharge at which complete 
lateral mixing takes place has been estimated to be 42.66 km. Transversi 
concentration distributions of tracers for the study segment of the Ottawa River 
are not readily available to check this estimated distance. However, it is 
reasonable to assume that lateral concentration gradients will persist for a 
considerable distance, although the presence of meanders and other turbulent 
conditions may result in the distance being less than 42,66 km. Thus, there exists 
a possibility for an overlapping of the convective period and lateral mixing zones. 
Unfortunately, analyses for this situation, either theoretically or through field 
studies, have not been reported in the literaturi. Fortunately, in the Ottawi 
River study segment, the convective period zone is not likely to have a significant 
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Figure 4. Tritium Distribution at the Deep River Water Supply Intake 



effect on the application of convectlon-drf fusion models, as per the analyses 
presented 1n the previous section. 

It Is recommended to further evaluate this aspect with the aid of 2-D model and the 
data on tracer distributions gathered during the September 1984 field studies 1n 
the Ottawa River study segment. If the convectlve period Is dominant within the 
lateral lixlng zone, then the measured concentration distributions would exhibit 
marked skewness. Based on such an evaluation, the need for further analyses could 
be deterii ned . 



3.3 Linikage Between Near- Field and Far-Field M^odlels 

3.3,1 IMear-Field and Far-Field Regions: 

In water quality modelling studies related to mixing zones. It is a common practice 
to consider the area of impact of receiving waters as near-field and far-field 
regions. Each of these regions 1s characterized by a distinct mixing process In 
the receiving water. Figure 5 shows a schematic representation of the various 
zones within these two regions below a bank outfall In a river channel, as well as 
the type of model applicable to each zone. 

The near-field region, closest to actual outfall, is predominantly influenced by 
the outfall configuration, jet-effect (I.e. momentum of the discharge), as well as 
bouyancy effects caused by density differences between the effluent discharge and 
receiving water. 

The far-field region, located some distance away from the outfall. Is characterized 
by passive convection and diffusion mechanisms due to ambient flow regime, as well 
as decay of pollutants, yithin the fir-field region, concentration gradients may 
be present In the vertical, literal and longitudinal directions, the extent of each 
being dependent on the convection, diffusion and other hydraulic processes in each 
direction. Generally, the longitudinal distances required to achieve complete 
vertical and literal mixing are in the range 50-100 times depth and 50-100 times 
width of channel, respectively! for the Ottawa River segment below the RoTphton 
NPD discharge, the corresponding distances are estimated to be 0.345 and 42.66 kii 
(Gore Si Storrie Ltd., 1983). Since the zone of vertical mixing in river channels 
tends to be relatively smiall (in comparison to the lateral mixing zone), the 
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linking of near-f1i1d and 3-D models 1s not conslderid In this study. 

3.3.2 Linkage Between IMear-Field and Far-Field 2-D Models: 

In order to develop a procedure for Unking the models appllciblt to near-field and 
2-D transverse rnlxing zone models, the near-field model developed by Ditmars (1969) 
and the far-field 2-D models, MIXCALBN and NIXCADIF, have been considered as 
examples. The theoretical aspects including mathematical formulations of the 
near-field «odel can be found In various publications (Hamdy, 1981; Ditmars, 
1969). The far-field models are based on the mathematical formulations presented 
in varioys publications (Gowda, 1980; Yotsukura and Cobb, 1972); these models are 
operational on the MOE computer system. 

From the predictions of the near-field model, concentration distributions at 
various distances downstream of the outfall can be obtained. These concentration 
profiles generally resemble the Gaussian distribution. The nearest distance from 
the outfall (xn) at which the near-field model predictions exhibit vertical 
uniformity can be determined from these predictions. 

The far-field 2-D models can then be utilized to predict the concentration 
distributions at distances equal to or greater than x^. For the far-field 
models, the nondlmensional transverse diffusion factor (6) applicable to the 
distance {n^) should be determined by a model calibration procedure so that the 
lateral concentration distributions at the distance (xm) predicted by the 
near-field and far-field models are in reasonable agreement. This empirical 
procedure will ensure a continyity 1n predicting concentration distributions in the 
near-field and far-field zones. 

3.3.3 Linkage Between Far=Field 2-D and 1-D Models 

The far-field 2-D model is applicable to the zone within which transverse 
concentration gradients are predominant. Thus, when the lateral concentration 
disbribytions approach yniformity, the 1-D far-field model can be applied. 
Theoretically, it is possible to estimate the longitudinal extent of the transverse 
mixing zone below an effluent discharge (Gowda, 1980 and 1984); for example, the 
theoretical distance in the Ottawa River channel below the Rolphton NPD discharge 
1s estimated to be 42.66 km, as stated earlier. However, because of the presence 
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of leandirs, etc., the distance could be smaller than 42.66 km. Therefore, 1t is 
bitter to chick the predicted lateral concentration distributions at various 
distances for the attalnaent of unlforiiity, and to determine the nearest distance 
(xj) below the outfall at which the uniformity of concentrations is achlivid. 
From a practical perspective, the distribytions can be considered to be nearly 
uniform when the absolute percentage difference between the minimum and maximum 
concentration levels reaches an arbitrary limit of 10 to 2011 as syggested by Holley 
et al (1972). Such a procedure does not significantly affect the predictions and 
will resylt in lessir computatlonai time requirements. 

The procedure for linking the far-fiild 2-0 and 1-D mode Is outlined above can be 
easily incorporatid for Instantaneous, flnite-timi as well as continuous discharge 
cases. 
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W SUMMARY AND CONCLUSIOIMS 

As part of the Ottiwa River nuclear spm contilngency model development studies, an 
analytical solytion for the 2-D f1mlte-t1me case has been derived. The derivation 
is based on the stream tube concept wherein the lateral co-ordinate Is represented 
by the cumulative discharge. An approximate Integration procedyre, known as the 
Laplace method, has been utilized in the derivations. The antlytical solution 
includes the effects of side-wall reflection of the material from channel banks. 

Expressions have been obtained for the expected arrival time and duration of 
passage of spilled material cloyd at a given station distant (x) from the source. 
Expressions are also presented for computing moving average values of 
reach-dependent parameters which include channel widths, depths and velocities, 
longitudinal dispersion coefficients, lateral diffusion factors and pollutant decay 
rates. 

A computer pogram for the 2-D finite-time model, TWODIFIN, has been developed in 
FORTRAN language. The documentation for the program presented in Appendix B 
Includes the program listing, input/output variables and a sample input and ouput. 

The importance of the convective period in the Ottawa River below the Rolphton NPD 
discharge has been assessed with the aid of the tritium concentrat1on-versus-t1me 
distributions monitored at Deep River, Petawawa and Pembroke during August ^ 
September 1981, following the tritium spill from the Rolphton NPD in August 1981, 
The convective period is found to dominate for about 16-20 km below Rolphton as 
indicated by a skewed concentration distribution at Deep River; however, the 
degree of skewness relative to that for a normal distribution (skewness = 0) is 
found to be small. Thus, It 1s concluded that the effect of convective period can 
be neglected 1n the modelling studies. It is recommended that this aspect be 
further examined with the aid of data from the field studies conducted during 
September 1984 and modelling studies. 

The relation between the convective period and lateral mixing zones has been 
examined. A detailed assessment could not be made because of a lack of relevant 
Information (either theoretical or field study data); it Is recommended to further 
evaluate this aspect with the aid of 2=0 model and the field data gathered dyring 
the September 1984 field studies. 
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Methods for linking near=f1eld and far-field models have been presented, A 
schematic representation shows the various zones 1n the near-field and far-field 
regions, and the models applicable to these zones. The methods for linking of 
iodels presented hirein will permit a continuity in applying the models for the 
near-field and far-field regions. 
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A.1 Evaluation of Integral by the Laplace Method 



The integral appearing in Equation 13 is 

ti 
I'f expT-f e \ - n (t - t')1 dt' 

ti 



/ ' ^^^ ( - tf - ^^^) ^ 



which can be rewritten in the form 



exp LA. nil J J 

T 



(A.l) 



Let (t - t')/t - T (A. 2) 

Thin, df - -t.dT; t - t' - tT (A. 3) 

The integration Units wiTI be as follows: 

As V .-ti, Ti = (t ' ti)/t (A. 4a) 

As V ►tj. Tj = (t - tj)/t (A, 4b) 

With this change of variables. Equation A.l becomes 



(A, 5) 



I - (-1) f " exp [X.li(T)] dT (A. 6) 



where X « c/t ; (A. 7a) 

h:(T) ^ =/l + T \ (A. 7b) 

a' = U/nt^) (A. 7c) 

M 



For larfe values of X, the major contribution to the viloe of the integral I 
arises from the Immediate vicinity of those points of the integration intirval 
(Ti, Tij at which h(T) has the largest values. The points at which h(T) 
attains maxima are obtained by dlffirentiating Equation (A. 7b) and setting the 
first derlvitivi to lero as follows: 

h'(T) = d [h(T)] - (A. 8a) 

IT 

- ~/-l + 1 \= (A.8b) 

Thust wi get T ^ a (A. 9) 



For h(T) to attain a maximum value, the second derivativi of Equation (A. 7b) should 
be negative. This, in fact, is the case as shown bellow: 

h"(T) - d^' [h(T)] - -2 (A.IO) 



Because of the fact that the major contribution to the value of the Integral, I, 1s 
in the vicinity of T = a. Equation A. 6 can be written 1n the approximate form 



a + e 
I = (>i) exp [X.h(T)] dT (A. 11) 

•'a - e 

in which e 1s some small positive number. 

Now, let us define 

s^ - h(a) - h(T) CA.12a) 

Then dT = ~2s.ds (A. 12b) 



IS 



The integration limits are as fonows: 



Lower H«1t: Sii = [h(a) - h(a-E)]^/^ (A. 13a) 

Upper limit: sj = [h(a) - h(a+e)]^/^ (A. 13b) 



Note that both S| and Sj are positivt since h(a) is the maximum valus 

of h(T). Substitution of Equations A. 12 and A. 13 in Equation A. 11 results 1n the 

fonowlng expression for the integril: 



I ixp [Xfh(a) - sm 2.s.ds (A.14) 

I T.h'(T) 



From Equations A. 9 and A* 10, 
Lt fl\ ^ 1 



i^lfiTTj) 



- [-2h"(a)]^l'/^'' - -Va^/2 



From Equations A, 14 and A, 15, we get 



,.h(a)] 



(A.lSa) 



(A. 15b) 



I = -VcT exp [X.h(a)] | exp (-Xs^) ds (A. 16) 

Si 



From Equations A. 7b and A, 9, 

h(a) - - (2/a) (A. 17) 



S# 



Now, 1st us ixpress ths 1ntigrat1on limits s^ and sj in ternis of T| and 
Tj. Substitution of Equations A. 7b and A, 17 "in Equation A, 12a results In 



^ = [h(ci) - h(T)] .= /-2 + 1 + T \ - (T - g) ^ 



(A. 18) 



Then, Si| = Tj - ot 



oTTf 



(A. 19 a) 
(A. lib) 



The integral appiarlng on the right hand side of Equation A. 16 can now be evaluated 
as foTlows 



/ 



exp (-Xs^).ds - 1 fF [erf(svT')] 



M^ I erf ( s j ^) . erf ( s .. /%)] 



Si 



w I erf I 



iMJhl- '"h^)M 



(A.20) 



Sybstituting the expressions from Equations A,17 and A.20 in Equation A,16, the 

following expression for the Integral, I, is obtained: 



I ^ - h[a#exp (-2X/ 
^l4X 



ii(^)j?l-i(Mj^ij'"- 



21) 



Frow Equations A. 4 and A. 7, 



erf 



[(^)i^i--iJ^-^) 



(A. 22) 
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Simnar ixprssslons for the other error function term of Equation A. 21 can be 
obtained. Thirefore, the final expression for the integral appearing In Equition 
13 bi comes: 

I = l? Bxv(-2^{Zn) ferf / /T - ./nJTj - ^'"^(jV'- ^)] (A.23) 



A.2 Derivation of Expression for Time of Paiiage 

Let us consider Equation 8 with K^j - 0, and define a relative concentration as 

^' = Cl/Cpgak - ci (2TftvfxE^/Wi u) (A. 24) 



1n which Cpgj|^ = peak concintratlon which occurs along q = q^ at x - ut. 
As stated in Section 2.1.4, Equation A. 24 gives an approximate ratio of a disired 
mlnirayip concentration to the peak concentration at a given location. From 
Equations 8 and A. 24, we get the expression 

f ' = expf~ (x - ut) ^ - JQ,- q) ^| (A,25) 

I 4Ext 4Eqt J 



Taking logarithms to base i and simplifying, we get the following quadratic 
equation 1n t ; 

t^ - [(2x/u) = (4 Ex/U^) In f ]t + {(x/u)^ + (E^/E^) [(q^ - q)/ul^l = (A. 26) 



For simplicity 9 we consider ix = 1 and my = 1, and substitute 

Eq ^ (BQ^u/b) from Equation 15a. Then, the solution of Equation A.26 is given 

by 

Equation A. 27 Is the same as Equation 17 given in Chapter II. 
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IB.1 Description of Variable Names 



TITLE 

m 

if 
IT 

Rqs 

THETA 

D'LT 
PS IP 

fl 

m 
i 

EXD 

BETA 

BP 
ZP 
UP 

AS(I.l) 

AS(I,2) 
AS(I,3) 
AS(I.4) 
AS (I, 5) 



Title of study. 

Number of transects (including the first transect at the source) 

Number of elemental strips at transect 

Number of lateral data points = NQ + 1 

Number of Input time- intervals 

Cumulative discharge between reference bank and pipe outfall 
location In the cross-section expressed as a friction of 
river discharge below the outfall. 

Temperatore correction factor associated with the decay 
rate (K). 

Value of 6 for estimating the number of images. 

Value of f ' for estimating passage times of cloud. 

Input times from source or outfall (IT values), hours. 

Input discharges from the outfall (IT values), raVs. 

Input concentrations (IT valyes). 

Distance from outfall to a given transect (NX valyes 
including a transect at X - 0), meters. 

Longltydlnal dispersion coefficients (NX values) for 
successive reaches of lengths (xj ~ xj_|), m^/s 

Non-dimensional transverse diffusion factor values for the 
successive reaches (NX values). 

Width exponent of the Leopold-Haddock equation for each reach. 

Depth exponent of the Leopold-Maddock equation for each reach. 

Velocity exponent of the Leopold-Naddock equation for each 
reach . 

Riference temperatyres (°C) 1n successive reaches at which 
the decay rates are known. 

Reference streamflow rates in successive reaches, mVs 

Reference channel widths in successive reaches, meters. 

Reference channel depths In syccessive reaches, meters. 

Reference channel velocities In successive reaches, m/s. 
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AS(I,i) Reference decay rates (base e) 1n successive reaches, per second. 

AO(Itl) Deilgfi temperatures (°C) in successivt reaches. 

AD(I,2) Design streamflow rates 1n successive reaches, m^7's 

AD(I,3) Channil width at the design streamflow, meters. 

AD (I, 4) Channel depth at the design streamflow, streamflow, 

AD(I,5) Channel velocity at the design streamflow, m/s. 

AD (I, 6) Decay rate coefficients at the design temptratures, per second. 

i Moving average value of channel width for a channel length 

Xi (from outfall to transect, i), meters. 

I Moving average value of channel depth, meters. 

i Moving average value of channel velocity, m/s. 

■ii Moving average value of decay rate coefficient, per second. 

EXX Moving ave. value of longitudinal dispersion coefficient, m^/s. 

if Moving average value of the transverse diffusion factor, 

TMl Tirae of arrival of cloud (ta), hours, 

TM2 Oeparture time of cloud {%)» hours. 

TMR Duration of cloud passage (tp), hours. 

■^ff Partial cumulative discharge (q), m^/%, 

IT River discharge, w'Vs 

m Ratio (QY/QT). 

if Mass flux input during a specified time-interval. 

WIN Totil mass flux input from the source* 

Tli Time of concentration prediction (t), hours and days. 

61 Concentration at a point (x,q) at a time, t. 

Vise Kinematic viscosity of water at the reference temperature, m-/s, 

VIST Kinematic viscosity of water at the design temperature, m^/s. 

FLAGEX Flag for ex', 1 for direct iniput and 2 for calculating 
fix from BansaTs relationship. 
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1.2 List of TWODtFIN Program 

C TWOOIFIN: 2-0 FINITE TIME FAR-FIELD MODEL FOR VERTICAL LINE SOURCE IN RIVERS 

C PROGRAM OEVELOPID lY: DR. T. P, H. 90WDA* WATER RESOURCES OIVISIONi 
C eORE f* STORRIE LTD. , TORONTO, ONTARia, 

C FOR THE WATER RESOUCES iRANCH. 

f ONTARIO MINISTRY OF THE ENVIRaNMENT, TORONTO. 

i' SEPTEMIER 19S4 

C « » ♦ » # ♦ » * * MAIN PRO© R A M # ♦#«##» # » 

6 

OIMENSiaN TITLEC20K, CTC6, 21, SI ». NTMC6K IM1C21 I, IM2(21 }. KTM(21 ) 
REAL*1 EKK C 6 K, A I , A2, A3 , A4 , AS , A4 , AD C 6 , 7 1 , AS C 6 , 7 1 , SUMT , WT C 1 OJ , 
t;S.lETA(6) , RQCgl > . K(6), RiK(A> . iP, IP, UP, aRTO, QSO, P^i RGB, TINi TMR,, 
21 X 1 1 TM 1 , TM2, AP 1(21 , 30 > ,• AP2 C 2 1 , 30 ) , AP3 ( 2 1 , 30 ) , AP 4 ( 2 1 * 30 ) , TOT , 
31 1 C 21 , 30 } . i2< 21 , 30 1 , 13 C 21 , 30 > r 14 C 21 , 30 ^ , Dl , D2« D3, D4, SMEXD C 6 I . 
4C1* C2, C3, C4, BEQ, TI» TJ, lETAC4)i TMC6, SI > , TECIOKp E)CDC6? , AMID* DLT, 
SISUM •( 4 ) , SUMK, ( 6 )• , STOT C 6 ) , VOL i^h, THETA,. STIC 51 , 30 J ,^ST J < 5 1 , 30 > ,• X X , 
6PRN (21 K PRP C21 ) , QE C 10 ) , CE ( 10 } ,. PSP (21 ) ,, PTP (21 ) , REYNQ., TR ISr TFAL,. 
7TR 1 .. TR2,, DP VL, PS I P .• V I SG , V I ST. i T ( 4 1^ i i ( A t , I ( 4 ^r U ( 6 } 
C 
Q DATA INPUT 

iREAOd, lOJTITLE 

READ( 1, *} NX. NQi IT, ROS, THETA. DLT, PSIP. VISC. IFLAG 
DO 20 Is I, IT 
20 READ(lr«)TE(I),iE(ir, CECII 
C READ HYDRAULIC, DECAY RATE AND TEMP. DATA FOR EACH TRANSECT FOR 
C REFERENCE AND DESIGN CONDI TiaNS 
00 22 1^1, NX 
SO TOCOO, 321, IFLAS 
30 READ C I , « ) X ( I) , EXD ( r ^^ BE'TA ( I ) ,, iP , ZP , UP 

30 TO 34 
32 READ 1 1 , « I X ( I K lETA ( U , BP , ZP , UP 
34 READ CI I * > ( AS C I , IS ) , IS^^I , 6 } 
READ C 1 ., « > C AD ( I , IS J , IS-1 , 2 ) 
C SCALE-UP HYDRAULIC PARAMETERS AND DECAY RATES TO DESIGN CONDITION 
PNO=FLQATCNQ) 
PQ-BILECI. O/FNOJ 
mm TO-AD ( 1 1 2 > / AS ( 1 , 2 ) 
AD (1,3 } =AS C 1 , 3 > •ORTO*»BP 
AD (1,4 ) =AS ( 1,4^ •QRTQ«*ZP 
AD ( 1 , 5 ) =AS C I , S J »@RTO*#UP 

AD C I ,^ 6 > ^^AS ( 1 , 6J •THETA** ( AD C 1 1 1 i - AS ( 1 , 1 H 
QQ T0C22, 36J, IFLAG 
C CALCULATE LONO 'L DISP 'N COEFF. USING iANSAL'S EON. 
34 DP VL=AD ( 1 , 4 > #AD C I , S } 

V I ST=V I SC » 1 . 029** C AS C I i 1 J - AD ( 1 , 1 H 
REYNO-DPVL/VIST 

EXDCIMC70SOO00. 0i#DPVL/REYN0*#O. 743 
22 CONTINUE 
C CALCULATE INPUT MASS FLUX FDR EACH TIME STEP AND TOTAL INPUT FLUX 
IDT^IT-1 
HIN*0. O 
00 24 L=l, IDT 
LL=L+l 

T IN^3600. 0* C TE ( LL ) -TE C L ) ) 
WT C L } =0. S# C OE C L ) *CE C L ) +C1E ( LL » *CE C LL ) ) 
24 WIN^WINH»-WTCL)»TIN 
NY=N0+1 
00 38 J=i.NY 



m 



PRNCJMSQS-RaCJ) 
30 PRP(JI-iQS+Ra(J> 

C COMPUTE MOVING AVERAQC VALUES OF HYBRAULIC PARAMCTERS. LONG 'L AND 
C LAT'L DISPERSION COEFFTa. . AHB BECA¥ RATE 

SUMK(l)sO. 

VOLCll'O. Q 

iSUMCll^O. 

STOTCll^O. 

SMiEXD(ll=0^ Q 

SBETACl*-0. 0' 

DO 76 I=3iNX 

11:^^1-1 

XX=X<IJ-X{I1) 

iX I'X X» { AD (1,3 1 *AD ( 1 1 „ 3 n 

iSUM( I ) =iSUM C I IJ +0. 5»i X I 

VOL C I J ^VOL ( r 1 > +0. 25«iX I* ( AD {1,4) +AD C 1 1 , 4 > > 

TOT*XX/AD(I, 51 

SUMK ( I ) =SUMK < 1 1 > + C TOT» AD C 1 , 6 H 

STOT ( 1 1 -STOT C 1 1 ) +TDT 

SMEXD C 1 1 =SMEXD C 1 1 1 •••KX»EXO C I J 

SIETAC I } -SBETAC 1 1 :^ + XX»iETA( I ) 

iCI)-iBUM(I)/XCI) 

Z ( I ) -VOL C I ) /iSUM ( i I 

UC I )-A0( l,r 2 J /{B C I l«ZC 1 1 J 

RK C I ) -SUMK C I ) /STOT ( I } 

EXXCI ) -SMEXBC II /X C I > 

BT(I)=SiETA(II/X(I) 

BE«1-B CI)/ ( IT C 1 1 »U C 1 1 > 
C CALCULATE VALUES OF TERMS DEPENDENT ON XCM 

A 1 « (1.0 / { 7. QB9Q 1 :i4#AD (1,211) »DSQR T ( C B ( I ) ♦U C I ) ) / ( EX X C II »B T ( I ) ) ) 

A2-#. 0*EXX<I) 

A3sDSQR T C ( U C I ) «U C I > / A2 ) +RK ( I ) ) 

A4-X(I)#*2/A2 

AS-U ( I ) «X < I ) / C 2. 0#EXX ( I ) ) 

A&-Q, 25* C A5+DLG# ( 1 . 0/DLT ) ) ##2/ ( A3»A3 ) -A4 

IF(A6. LT. O. 01 A6-0. O 

AMiD^DSeRT ( A6/BEQ ) 
C INITIALIZE liMAGE TERM ARRAY VALUES TO ZERO 

00 70 J=t, NY 

IMl(J)-0 

IM2(J)-0 

DO 70 N=l,30 

APKU, N)=0. 

AP2(J, N)=0. D 

AP3CJ, N)*0 

AP4(J. N)»0. 

ilCJ, N)=Q. 

i2CJ,N)*0. 

iSCJ, N)=0. O 
70 i4(J,N)-0„0 
C COMPU'TE AND STORE VALUES OF IMASE TERMS liNDEPBMBENT OF TIME 

DO 74 J=1,NY 

C ALL R I MASE C AMi D , PRN , J, M 1 , M4 ) 

C ALL R I MA§E C AMiD,, iPRiP , Ji M2 , M3 ) 

IMi(J)=AMAX0CMi,M2l 

IM2{J)-AI1AX0(M3. H4) 

MJI-IMICJ) 

Hy3'=If12CJ> 

PSPCJI--PRiN(JI 
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PTPCJ)--PRPCJ) 

DO' 72 N-1. HJl 

CALL REFLIC CPRN, BEQ, A3* A4i AS, J. Nr API, 81 } 

CALL REFLECCPRP. iEQ» A3,r A4.» A5, Ji N. AP2, B'2) 
73 CONTINUE 

BO 74 N«2, M^ 

CALL RiEFLECCPTP, WEQ. A3,, A#. AS, J, N; AP3, 13) 

CALL REFLEC CPSPr lEQi A3i A#, AS, J, N, AP4, B4 ) 
?4 CONTINUE 
GENERATE PASSAGE TIMES TMCI,K^ AND COMPUTE VALUES OF TIME-DEPENDENT TERMS 

CALL PASTIMCX, U, EXJCi PSIP. I, TMI, TMg) 

TMR^TM2-TM1+TEC ITI-TE( 1 } 

IFCTMR. LE 10. OJ NT^20 

IFCTMR. QT. 10. 0. ANC. TMR. LE. 40. O) NT=40 

IFCTMR. OT. 40. 0) NT^SO 

TRIS=CO., e#X( I I/C3600. 0*U< I ) ) l-TMl 

TFAL^TMR-TRIS 

TFN^FLOATCNT) 

TFN-sO. 4#TFN'i-0. S 

NT1-^IFIX<TFN) 

NT 1^0. 4#NT 

NT2-SNT-NT1 

TRl^TRIS/CNTl-lJ 

TR2=TFAL/CNT2-li 

00 26 K-1,NT 

IFCK. LE. NT! } TMC I, KJ=TH1+TR1*CK-1 1 

IFCK. QT. NTl I TMC I, KI^TMCI, NTl H-TR2«CR-NT1 I 

DO 26 L-1, IDT 

LL=L+1 

T I ^3600 . 0* ( TM C I , K I -TE C L > I 

T J=3600 . 0« C TM C I , K » -TE ( LL I ) 

IFCTI. LT. 0. 0. OR. TJ. LT. 0. 0)G0 TO 40 

STICR., D^OSQRTCTI) 

STJCK, L)=OSQRTCTJ^ 

SO TO 26 
40 STKK. L)=-99. 9 

STJCK, L)=-99. 9 
26 CONTINUE 
CALCULATE CONCN'S FOR EACH POINT (I, J) AT EACH TIME TMCI,K) 

,B0 38 J=1,NY 

KTMCJI^O 

Ml-IMICJ) 

M2=IM2CJ) 

DO 96 K-1,NT 

SUMT^O. 

DO 78 L*lf IDT 

IF(STI<K, L). LT. 0. Q. OR. STJCK, L^. LT. O. OIQO TO 84 

C1=0. O 

G2«0. Q 

C3«0. Q 

C4«0. O 

DO 10 N=i,. Ml 

CALL FUNER (API, STI, STJ, A3. Jt N, K. L, Dl > 

CALL FUNER CAPS, STI. ITJ, A3* J. N, K, L, D2) 

C1=C1+B1CJ,NJ»D1 
80 C2^C2+i2 C J, N I *D2 

00 S2 N*2, M2 

CALL FUNER C AP3, STI, ST J, A3, J, N, K, L, D3 > 

CALL FUNER C AP4, STI. STJ, A3, J. H, K. L, D4) 

C3sC3+e3CJ, N)»P3 
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,82 G 4^C 4+1 4 ( J . N ) # 54 

lUMT-SUflT+yT ( L ) ♦ ( C 1 +C 2+C3+C 4 ) 
7e CONTINUE 
S4 C T ( I / J. K » -A 1 «SUMT 

IFCR. LE. 19. QR, K. EQ, NTiaO TO 86 

IFCCTd. J. K J. @E. 0.000011 00 TO 86 

KTMCJ^^K 

00 TO isa 

86 CONTINUE 

KTM ( J > ^NT 

eS CONTINUE 

NTMCI>--1 

00 28 J=l, NY 

IFCNTM( r > . GT. KTM( J) > SO TO 2B 

NT«CI^*KTI1CJ) 
28 CONTINUE 

NT-NTH CI)' 

PO 76 J=l. NY 

KM-KTMCJ^+l 

IF (KM. LE. 20 J QQ TO 76 

IFCKH. 3E. NT^ GO TO 76 

DO 42 RsKMi NT 
42 CTCI* U,r K)-0^ 
76 CONTINUE 
C P^'INT CONC'N DISTR'N AT X(I,JJ AT THE TIMES TM(I,R> 

iyRITECI, 10) TITLE 

lYsNY/lO 

IFCIY. LE. II lY^l 

DO 44 1=2, NX 

MT=NTM(0 

ISTN^'I-1 

MRITECi.52> ISTN 

WRITEC1.12) 5C( I), A:DCI. IK ADCI, 2)^ 

W^R ITEC 1 ,141 e < II , Z ( I ) . UC I ) , EXX ( U , ITC I } r RKC I ) 

Un ITE { I * 16 I CRQi C J } , J=l .' NY, I Y ) 

DO 44 K-1, NT 

TDAY=TM(I,K»/24,. O 
44 WR ITE ( 1 . li ) ¥\. TM CI , R ) , TDAY. C CT ( 1 1 J, K. > . J- 1 . NY. I Y I 
C FORMAT STATEMENTS 
10 FORMAT C20A4> 
12 F0RMATC3X. 'X-S FlO. I, 3Xi ^TEMF»', F6. 2, 3X, 'Q^% FS. 2) 

14 FORMATOX. "'AW,. VALUES: '. 2M. 'B= ', F7., 2. 2X. 'Z- % F6. 2, 2X, 'U« ',, F6. 3. 
«3X. 'EX^', FS. 3, 2X, ^iETA^% F7. 4. 2X. "^D-S ElO. 4/) 

16 F0RMATC2X, 'NO. '. 1X« 'TIMECHRS. : DAYS ) % IKlX.Ff. 2)/) 

15 FORMATCIX, 13. 1X,_F8., 2, 2X.^F5,. 1, IX, llClX,Ff. 4) ) 
52 F0RMAT(/3X. 'STATION ',121 

STOP 
END 

C »##»*# S U 1 R y T I N E S #**### 

C 

C CALCULATE ERROR FUNCTION VALUES 



C 



SUBROUTINE FUNER CAP« SI, SJ, AR, 14, M, Kl, LI , D J 

ilAL»8 AP ( 2 1 . 30 1 , S I C i 1 , 30 ri S J C 5 1 , 30 ) . AR . D. 1 1 , i J 

IF ( S J CICl , L I) ., EQ . Q. Q I ©0 TO 1 

iU^OERFC CAPCK, M) /SJCR1« Li ) :i-AR»SJCKl, LI ) ) 

QQ TO 12 

10 iJ=i.o 

12 IFCSICK1,LII. EQ. 0. 0)GO TO 14 
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ai^OERFC CAPCI^. MK/SI CKI. LI I l-AR«SI CKl, LI ) I 
©0 TO 16 
14 11=1.0 
16 O-lJ-il 
RETURN 
ENO 
€ 
C EST I MAT I ON OF PASSAGE TIMeS OF CLOUS 

C 

SUlROUTIiNl PASTIMCXI, UI, EXI, PSI» IK. TXl, TX2) 

REAL#e X I C 4 1 , UI ( 6 ) . EX I ( 6 ) . TX I , TX2, PSX, UX . UT. XUE. SXU, PS I 

PiX-EX I ( I X ) «0L08 C 1 . 0/PS I ) 

UX-UI(IXI»XI(1XI 

UTsU I C I X 1^ ♦U I ( I X > #3600 „ 

XUE-UX+2. 0#PSX 

sxu^psx*cyx+psx) 
sxu*2. o«DeaRT(sxuj 

TXl^CXUE-SXyi/UT 
rFCTXl.LE. 0.01 TXl^O. 05 
TX2*0. S# ( XUE+eiU I /UT 
RITURN 
END 

c 

€ COMPUTE VALUES OF IMAQE-DEPEMDENT TERMS 

€ 

SUIRQUTINE REFLEC (PR* BEi ETA. ALPA, AX„ M. K, ANi BH) 

REAL«e PR iZl > . ANCaii 30> . BNCgl, 30) i AR. iEi ETAi ALPA, AX, PN 

PH^ i 2,. 0« C K- 1 ) +PR C M H ♦»2 

AN C 11 1 K ) ^BS@RT ( ALPA+0. 2S#iE»PN ) 

AR=ETA*ANCM, K» 

iN c Ml K ) =DEXP ( AX-2. 0*AR ) /OSQRT ( AR I 

RETURN 

END 
C 
C CALCULATE THE NUMllR OF IMAOES 

g- 

SUBROUT I NE R I HASE C AM . PR . J i MT t , MT2 J 

REAL«8 AM, PR 'C 21^ 

ANl ^-0. 5»PR C J ) ~AM-0. S 

AN2^-'0. 5*PR C J } +AM+0. 5 

Nl^INTIANl) 

N2^INT<ANgl 

MT I = I AiS C N II + 1 AiS C N2 Y 

AN3^0. 5«PR(J)-AM-0. 5 

AN4-0 . 5«PR < U ) + A«+0. S 

NI3=INT<AN3) 

N#-I!NT(AN4) 

MT2*IAiSCN3)+IAiSCN4) 

IFCMTl. LE., 1* MTl = l 

IFCMT3,. LE. 21 MT2^2 

IFCMTl. ©T. 301 MTli^OO 

IF(MT2. ©T. 30) MT2^30 
RETURN 

END 



B.3 input Data Sequence 



The sequencs of data Inputs ari as follows; 



Variable 



Format 



Remarks 



TITLE 

NX, NQ, IT, RQS, THETA, 
DLT, PSIP, FLAGEX 

TE(I), QE(I), CE(I) 

If (FLAGEX - 1): 

X(I). EXS(I), BETA(I). BP, ZP, UP 



If (FLAGEX - 2): 

X(I), BETA(I), BP. IP. UP 

AS(I, IS), IS = 1,6 

AO(I, ID), ID - 1,2 



20A4 


Oni lini 


Free 


One line 


Free 


IT number of lines 


Free v 


Each set of these 

three lines is 

intered NIX timis* 


Free f 




Free j 




Free 





The program is set up to read the Input data from an auxiliary device such 
as 1 disk or a tape. 



m 



e.4 ixatupie of Inpyt/Output 



Inpyf Data 

TRITIUII SPILL FROH ROLPHTON NPO TO OTTAWA R. : MJQ. 17-25, IfSl 
5 5 f 0. 1.03 0.001 0. 1 0. OOOOQIQIOI 2 

0. 10. o. 

1. O 10. 130. O 
24. 10. 740. 
48. 10. ?20. 

72. 10. eio. 

f6. 10. O §S0. 
120. O 10. 400. 
168. 10. SO. 
192. O 10. O O. 

O. a O. 001 O, 375 O. 375 0. 25 
20.0 2¥8. 7 1170.0 16.21 0.0157 0.0 
20. Q 493. 1 

15788.0 0.001 0.375 0.375 0.25 
20.0 298.7 1170.0 16.21 0.0157 0.0 
20. O 493. 1 

27295. 0. 001 0. 375 0. 375 0. 25 
20.0 298.7 1016.0 18.49 0.0159 0.0 
20. 493. 1 

45472. 0. 001 0. 375 0. 375 0. 25 
20.0 298.7 1800. 15.65 0.0106 0.0 
20. 493. 1 

62292. O 0. 001 0. 375 0. 375 O. 25 
20. O 298. 7 1777. O 4. 36 O, 0386 0. Q 
20. O 493. 1 



ii 



Oytpiit 

TRITIUn SPILL FROM ROLPHTON NPD TO OTTAWA R, 



AUQ. 17-23 1 1981 



STATION 1 
















X- 


1578S. 


TEMP* 


20. 00 Q^ 


493. 10 










AVO 


VALUES: 1 


-1411. ' 


?6 Z= 19. 56 


U=s 0. 018 


EX- 148. 


064 1ETA= 


0. 0010 


KD-0. OOOO 


NO. 


TIME(HRS. 


:DAYS> 


0. 00 


0. 20 


0. 40 


0. 60 


0. 80 


1. oo 


1 


36. 73 


I. 5 


i. 3749 


0. 0014 


0. 0000 


0. 0000 


0. OOOO 


0. OOOO 


m- 


45. 60 


1. 9 


2. 7719 


0. 0100 


0. 0000 


0. 0000 


0. OOOO 


0. OOOO 


& 


34. 48 


2. 3 


5. 4843 


0. 0379 


0. 0000 


0. 0000 


0. OOOO 


0. OOOO 


« 


63. 36 


2. 6 


S. 7444 


0. 1013 


0. OOOO 


0. OOOO 


0. OOOO 


0. OOOO 


i' 


72. 24 


3. 


12. 6790 


0, 2196 


0. 0000 


0. OOOO 


0. OOOO 


0. OOOO 


6 


81. 12 


3. 4 


17. 0205 


0. 4121 


0, 0000 


0. OOOO 


0. OOOO 


0. OOOO 


f 


89. 99 


3. 7 


21. 5409 


0. 6941 


0. 0002 


0. OOOO 


0. OOOO 


0. OOOO 


8 


98. 87 


4. 1 


26. 1783 


1. 0751 


0. 0003 


0. OOOO 


0. OOOO 


0. OOOO 


f 


107. 75 


4. 5 


30. 5418 


I. 5574 


0. 0011 


0. OOOO 


0. OOOO 


O. 0000 


10 


116.63 


4.9 


34. 4990 


2. 1366 


0. 0025 


0. OOOO 


0. OOOO 


0. OOOO 


11 


125. SO 


5.2 


38. 1958 


2. 8020 


0048 


0. OOOO 


0. OOOO 


0. OOOO 


12 


134. 38 


5. 6 


41. 2804- 


3. 5386 


0. 0088 


0. OOOO 


0. OOOO 


0.0000 


13 


143. 26 


6.0 


43. 7243 


4. 3283 


0. 0148 


0. OOOO 


0. OOOO 


0. OOOO 


14 


152. 14 


6. 3 


45. 2423 


3. 1521 


0, 0236 


0. OOOO 


0. OOOO 


0.0000 


13 


161. ©2 


6. 7 


43. 8284 


3. 9879 


0, 0357 


0. OOOO 


0. OOOO 


0. OOOO 


l& 


169. 89 


7. I 


47. 8750 


6. 3239 


0. 0319 


0. OOOO 


0. OOOO 


0. 000© 


17 


178. 77 


7. 4 


48. 3728 


7. 6337 


0. 0727 


0. 0001 


0. OOOO 


0. OOOO 


IS 


187. 65 


7. S 


43. 6028 


8- 4060 


0. 0994 


0. 0001 


0. OOOO 


0. OOOO 


if 


196. 33 


Q. 2 


48. 5121 


9. 1312 


0. 1296 


0. 0002 


0. OOOO 


0. OOOO 


20 


240. 20 


10. 


43. 7084 


1 1 . 3399 


0. 3636 


0. 0018 


0. OOOO 


0. OOOO 


21 


283. 86 


11. 8 


36. 7844 


12. 9605 


0. 7052 


0. 0078 


0. OOOO 


0. OOOO 


22 


327. 53 


13. 6 


30. 7593 


12. 3959 


I. 0863 


0. 0218 


0. OOOO 


0. OOOO 


23 


371. 20 


13. 5 


25. 8918 


12. 2139 


1, 4376 


0. 0463 


0. OOOO 


0. OOOO 


24 


414. 87 


17. 3 


21. 9793 


1 1 . 2732 


I. 7176 


0. 0915 


0. OOOO 


0. OOOO 


25 


458. 54 


19. 1 


18. 8062 


10. 2607 


1.9135 


0. 1243 


0, OOOO 


0. OOOO 


24 


502. 20 


20. 9 


16. 2027 


9. 2673 


2. 0305 


0. 1709 


0. OOOO 


0. OOOO 


27 


545. 87 


22. 7 


14. 0425 


a. 3346 


2. 0812 


0. 2175 


0. OOOO 


0. OOOO 


2B 


589. 54 


24. 6 


12. 2322 


7. 4787 


2. 0801 


0. 2608 


0. OOOO 


0. OOOO 


29 


633. 21 


26. 4 


10. 7019 


6. 7035 


2. 0404 


2986 


0. OOOO 


0. OOOO 


30 


676. 87 


28. 2 


9- 3983 


6. 0064 


I. 9733 


0. 3299 


OOOO 


0. OOOO 


31 


720. 54 


30.0 


8. 2805 


5. 3821 


1. 8877 


0. 3542 


0. OOOO 


0. OOOO 


32 


764. 21 


31. 8 


7. 3167 


4. 9244 


1. 7906 


0. 3715 


0. OOOO 


0. OOOO 


33 


807. 88 


33. 7 


6. 4814 


4. 3266 


I. 6871 


0. 3825 


0. OOOO 


0, OOOO 


34 


351. 55 


35. 5 


5. 7543 


3. 9826 


1. 5811 


0. 3977 


0. OOOO 


0. OOOO 


35 


895. 21 


37. 3 


5. 1190 


3. 4863 


1. 4754 


0. 3881 


0. OOOO 


0. OOOO 


36 


938. 88 


39. 1 


4. 3621 


3. 1326 


1. 3719 


0. 3942 


0. OOOO 


0. OOOO 


37 


982. 53 


40. 9 


4. 0723 


2. 8167 


1. 2719 


0. 3770 


0. OOOO 


0. OOOO 


38 


1026. 22 


42. 8 


3. 6404 


2. 3343 


1. 1765 


0. 3670 


0. OOOO 


0. OOOO 


39 


1069. 89 


44. 6 


3. 2586 


2. 2817 


1. 0861 


0. 3550 


0. OOOO 


0. OOOO 


40 


1113. 55 


46. 4 


2. 9204 


2. 0556 


1. 0010 


0. 3414 


0. OOOO 


0. OOOO 


41 


1157.22 


48. 2 


2. 6202 


1. 8330 


0. 9212 


0. 3267 


0. OOOO 


0. OOOO 


43 


1200. 89 


50. 


2. 3533 


1.6713 


0. 8469 


0. 31 13 


0. OOOO 


0. OOOO 


43 


1844. 56 


31. 9 


2. 1153 


1 . 5083 


0. 7777 


0. 2955 


0. OOOO 


0. OOOO 


44 


1288. 23 


33. 7 


I . 9034 


1.3619 


0, 7137 


0. 2796 


0. OOOO 


0. OOOO 


45 


1331. 89 


55. 3 


1. 7140 


1. 2304 


0. 6544 


0. 2637 


0. OOOO 


0. OOOO 


46 


1375. 56 


57. 3 


1. 5445 


1. 1121 


0. 5997 


0. 2482 


0. OOOO 


0. OOOO 


47 


1419.23 


59. 1 


1. 3928 


1 . 0056 


0. 5493 


0. 2330 


0. OOOO 


0. 0000 


48 


1462. 90 


61. 


1 . 2568 


0. 9097 


0. 5029 


0. 2183 


0. OOOO 


0. OOOO 


49 


1506. 57 


62. 8 


1. 1347 


0. 8233 


0. 4603 


0. 2042 


0. OOOO 


0. OOOO 


30 


1550. 23 


64.6 


I . 0251 


0. 7455 


0. 4211 


0. 1906 


0. OOOO 


0. OOOO 



m 



STATION 2 

Km 272^§ TEfiP* 20. 00 Q-= #^3. 10 

AVQ. VALUES : 1-1372. 7f 1= SO. 12 U= 0. Oie 



NO,, TtMECHRS. :DAVS» 



I 
2 
3 
4 
3 
6 
7 
8 
'f 

10 
11 
12 
13 
l# 
15 
16 
17 
18 

20 
21 
22 

23 
24 
29 
26 
27 
28 
29 
30 
91 
.32 
33 
34 

as 

36 

37 
38 

3f 
#0 
41 
42 
43 
44 
45 
46 
47 
48 
4f 
§0 



91 . 85 

105. 63 

119. 40 

133. 17 

146. 94 

160. 71 

174. 49 

ISO. 26 

202. 03 

215. 80 

229, 57 

243. 35 

257. 12 

270. 89 

284. 66 

298. 43 

312. 21 

325. 98 

339. 75 

387. 18 

434. 62 

482. 06 

529, #9 

576. 93 

624. 36 

671. 80 

719. 23 

766. 67 

814. 10 

861. 54 

908. 98 

956. 41 

1003. 85 

1051. 28 

1098. 72 

1146. 15 

1193. 59 

1241. 03 

1288. 46 

1335. 90 

1383. 33 

1430. 77 

1478. 20 

1525. 64 

1573. 07 

1620. SI 

1667. 9S 

1715. 38 

1762. 82 

11 10. 25 



7. 3 

8. 4 



O 

6 
1 
7 
3 
9 
4 

6 
2 
1 
1 
1 
1 

O 


9 



9. 

9. 

10. 

10. 

11. 
11. 

12. 
13. 

13. 
14. 
16. 
18. 

20. 
SS. 
24. 
26. 
28. 
30. 
31. 
33. 9 
35. 9 
37. 9 
39. 9 
41. 8 
43. 1 
45. 8 
47. B 
4f . 7 
51. 7 
53. 7 
55. 7 
57. 6 
39. 6 
61. 6 
63. 6 
65. 9 
67. 5 
69. 5 
71. 5 
73. 5 
75. 4 



0.00 

1. 3215 

2. 5017 

4. 1161 
6. 1004 
B. 3462 

10, 7215 

13. 1054 

15. 3675 
17, 4338 

19. 2618 

20. 8171 

22. 0750 
23. 0302 

23. 6978 

24. 1078 
24. 2966 
24. 3015 
24. 1568 
23, 8929 
22. 3945 
20. 4619 
13. 4478 

16. 5136 

14. 7266 
13. 1084 
11. 6592 
10. 3691 

9. 2244 
i. 2103 
7, 3124 
6. 5174 

5. 8132 
1190 
6352 
1436 
7067 
3182 
9723 
6640 
3892 
1439 
924© 
7290 
5538 
3971 
2567 

1. 1308 
I. 0180 
0. 9168 
0. 8259 



11, 

11, 



0. 20 

0. 0793 
O. 2022 
O. 4233 

0. 7669 

1 . 2:450 

1. 8549 

2. 5013 
3984 
2740 
1755 
0736 
9425 
7599 

8. 3079 

9. 1739 
9, 7512 

10. 2379 

10. 6362 

10. 9310 

4914 

4132 

10. 9625 

10. 3097 

9. 5632 

8. 7894 

8. 0275 

7. 2997 

6. 6176 

5. 9865 

3. 4073 

4. 8791 
4. 3993 

3. 9649 
3. 5723 

2181 
8989 
6114 
3526 
1 197 
9101 
7216 
1. 5519 
I. 3993 
1, 2619 
I. 1383 
1 . 0270 
O. 9268 
0. 8365 
O. 7552 
0. 6820 



EX" 130.245 iETA= 0,0010 KD=0. OOOOE 00 
0.40 0. m 0.80 1.00 



O. 0000 
0. 0002 
O. 0009 
O. 0028 
O. 0072 
0. 0137 
0. 0307 
0. 0543 
O. 0889 
O. 1363 
O. 1977 
2737 
3642 
4683 
3849 
O. 7110 
0. 8436 

0. 9859 

1. 1299 

1. 6182 

2. 0457 
2. 3751 
2. 5990 
2. 7265 
2, 7738 
2. 7582 
2. 6956 
2. 5994 
2. 4806 
2. 3477 
2. 2072 
2. 0639 
I. 9215 
1. 7823 
1. 6483 
1. 5205 
1. 3996 
1. 2860 
1. 1799 
1. 0811 
O. 9894 
O. 9047 
O. 8263 
O. 7546 
O. 6885 
O. 6278 
O. 5723 
0. 5214 
0. 4749 
0. 4324 



0. 0000 
0. 0000 
0000 
0000 
0000 
0000 
0000 
0001 
O. 0002 
0. 0005 
O. 0010 
0018 
0031 
0050 
0077 
0.0114 
0. 0163 
O. 0224 
0. 0300 
0. 0679 
O. 1238 
0. 1935 
0. 2701 
0. 3461 
4181 
4800 
3306 
5692 
5960 
O. 6120 
O. 6185 
O. 6167 
0. 6082 
O. 5942 
O. 5759 
0. 3544 
0. 9306 
0. 5052 
O. 4790 
0. 4524 
0. 4298 
0. 3996 
O. 3740 
0. 3493 
0. 3253 
0. 3028 
O. 2812 
O. 2608 
O. 2415 
O. 2234 



0. 0000 
0. 0000 
0. 0000 

o. oooo 

0. 0000 
O. 0000 
0. 0000 
0000 

oooo 

0000 

oooo 
o. oooo 
o. oooo 
oooo 
oooo 
oooo 

0001 
0. 0001 
O. 0002 

0009 

0027 
0062 

0119 
0200 
0303 
0424 
O. 0593 
O. 0691 
O. 0825 
0. 0931 
O. 1066 
O. 1165 
O. 1249 
O. 1316 
O. 1366 
O. 1401 
O. 1420 
0. 1426 
O. 1420 
O. 1404 
0. 1379 
0. 1347 
O. 1308 
0. 1263 
0. 1218 
0. 1168 
O. 1117 
0. 1063 
0. 1012 
O. 0960 



0. OOOO 

o. oooo 
o. oooo 
o. oooo 
oooo 
oooo 
oooo 
oooo 
oooo 
oooo 
oooo 
oooo 
o. oooo 
o. oooo 
o. oooo 
o. oooo 
0. oooo 
0. oooo 

O. 000© 

0. oooo 
o. oooo 
o. oooo 

000© 

oooo 
oooo 
oooo 
oooo 
oooo 

©©00 

o. oooo 
o. oooo 

0. ©000 

o. ©oo© 
o. oooo 
o. oooo 
0, oooo 
o. oooo 
o. oooo 
o. oooo 
o. oooo 
0. oooo 

000© 

oooo 
oooo 
oooo 
oooo 
oooo 
oooo 

©000 
000© 



STATION 3 
















X« 


45472. 


TEMP= 20. 00 Q* 


493. 10 










AVQ. 


VALUES: 8 


-1903.26 


Z- 20. 34 


U* 0. 016 


eX= 143. 


641 BETA- 


0. 0010 


m^O. 0000 


UQ. 


TIMECH8S. 


:DAY8J 


0. 00 


0. 2© 


0. 40 


0. 60 


0. 80 


1. 00 


I 


222. 49 


9. 3 


1 . 3289 


0. 3341 


0. 0062 


0. 0000 


0. 0000 


0. 0000 


2 


244. 94 


10. 2 


2. 0943 


0. 9949 


0. 01 96 


0. 0000 


0. 0000 


0. 0000 


3 


267. 38 


11. 1 


3. 0193 


0. 9303 


0. 0333 


0. 0002 


0. 0000 


0. 0000 


4 


289. 83 


12. 1 


4. 0439 


I. 3947 


0. 0629 


0. 0004 


0, 0000 


0. 0000 


S 


312. 27 


13. 


9. 1291 


I. 9129 


0. 1077 


0, 001 1 


0. 0000 


0. 0000 


6 


334. 71 


13.9 


6. 2096 


2. 4S34 


0. 1698 


0. 0022 


0. 0000 


0. 0000 


7 


397. 16 


14. 9 


7. 2402 


3. 0326 


0. 2907 


0. 0043 


0. 0000 


0. 0000 


8 


379. 60 


19. 8 


3. 1992 


3. 6871 


0. 3500 


0. 0077 


o. oooo 


0. 0000 


f 


402. 09 


16. 8 


9. 0490 


4. 2770 


0. 4666 


0. 0127 


0. 0001 


0. 0000 


10 


424. 49 


17. 7 


9. 7901 


4. 3360 


9930 


0. 0197 


0. 0002 


0. 0000 


11 


446. 93 


18. 6 


10. 4192 


9. 3324 


0. 7413 


0. 0292 


0. 0003 


0. 0000 


12 


469. 38 


19. 6 


10. 9264 


3. 3187 


0. 3930 


0. 0412 


0. 0006 


0. oooo 


13 


491. 82 


20. 9 


1 1 . 3298 


6. 2303 


1.0497 


0. 0561 


0. 0010 


0. oooo 


14 


914. 27 


21. 4 


1 1 . 6339 


6. 3863 


1.2081 


0. 0739 


0. 0016 


oooo 


19 


936. 71 


22.: 4 


1 1 . 8468 


6. 8866 


1. 3691 


0. 0949 


0. 0024 


0. 0000 


16 


999. 16 


23. 3 


1 1 . 9794 


7. 1332 


1. 9191 


0. 1178 


O. 0O34 


0. 0001 


17 


981. 60 


24. 2 


12. 0409 


7. 3292 


1. 6690 


0. 1433 


0. 0048 


0. 0001 


IS 


604. 04 


29.2 


12. 0404 


7, 4781 


1. 8040 


0. 1714 


0. 0066 


0. 0002 


If 


626. 49 


26. 1 


1 1 . 9864 


7. 3837 


I. 9337 


0. 2011 


0. 0087 


0. 0002 


20 


689. 91 


28. 6 


11. 6439 


7. 6843 


2. 2244 


0.2949 


0. 0163 


0.0006 


21 


744. 92 


31. 


1 1 . 0993 


7. 9331 


2. 4380 


0. 3716 


0. 0270 


0. 0014 


22 


903. 94 


33. 9 


10. 4377 


7. 3400 


2. 9768 


0. 4949 


0. 0406 


0. 0028 


23 


862. 96 


39.9 


9. 7183 


7. 0029 


2. 6491 


0. 5302 


0. 0963 


0. 0048 


24 


921. 98 


38. 4 


a. 9819 


6. 6085 


2. 6694 


0. 3943 


0. 0735 


0. 0079 


23 


980. 60 


40. 9 


8. 2940 


6. 1837 


2. 6366 


0. 6462 


0. 0912 


0. 0111 


26 


1039. 62 


43. 3 


7. 9926 


5. 7479 


2. 9728 


0. 6930 


0. 1087 


0. 0193 


27 


1098. 64 


49. 8 


6. 8878 


5. 3149 


2. 4830 


0, 7113 


0. 1232 


0. 0201 


28 


1197. 66 


48. 2 


6. 2690 


4. 3939 


2. 3748 


0. 7260 


0. 1403 


0. 0293 


29 


1216. 67 


90. 7 


9- 6867 


4. 4908 


2. 2944 


0, 7306 


0. 1333 


0. 0308 


30 


1279. 69 


93.2 


3. 1334 


4. 1095 


2. 1268 


0. 7263 


0. 1647 


0. 0363 


31 


1334. 71 


99. 6 


4. 6639 


3. 7921 


1. 9999 


0. 7147 


0. 1736 


0. 0417 


32 


1393. 73 


98. 1 


4. 2166 


3. 4193 


1. 8643 


0. 6972 


0. 1804 


0. 0468 


33 


1492. 79 


60. 9 


3. 8090 


3. 1112 


1. 7337 


0. 6730 


0. 1390 


0. 0916 


34 


1911. 77 


63. 


3. 4384 


2. 9272 


1.6103 


0. 6492 


0. 1877 


0. 0998 


3S 


1970. 79 


69. 4 


3. 1023 


2. 3663 


1. 4898 


0. 6209 


0. 1886 


0. 0999 


36 


1629. 81 


67. 9 


2. 7979 


2. 3274 


1.3790 


0. 3909 


0. 1879 


0. 0626 


37 


1688. 82 


70. 4 


2. 9229 


2. 1091 


1. 2664 


0. 3398 


0. 1837 


0. 0651 


31 


1747. 84 


72. 9 


2. 2736 


1.9100 


1. 1642 


0. 9283 


0. 1824 


0. 0670 


39 


1806. 36 


79. 3 


2. 0488 


1 . 7288 


1. 0689 


0. 4969 


0. 1780 


0. 0683 


40 


1969. 88 


77. 7 


1 . 8460 


1. 3641 


0. 9793 


0. 4660 


0. 1729 


0. 0690 


41 


1924. 90 


80. 2 


1. 6631 


1. 4145 


0. 8964 


0. 4358 


0. 1670 


0. 0692 


42 


19Q3. 92 


32. 7 


1.4982 


1 . 2783 


0. 8197 


4065 


0. 1607 


0. 0699 


43 


2042. 94 


89. 1 


1.3493 


1. 1338 


O. 7488 


0. 3784 


0. 1540 


0. 0682 


44 


2101. 99 


87. 6 


1. 2196 


1. 0443 


0. 6834 


0. 3516 


0. 1471 


0. 0672 


49 


2160. 97 


90. 


I. 0990 


0. 9433 


0. 6233 


0. 3261 


0. 1400 


0. 0698 


46 


2219. 99 


92. 9 


0. 9863 


0. 8922 


0. 9680 


0. 3019 


0. 1329 


0. 0641 


47 


2279. 01 


99. 


0. 0889 


0. 7697 


0. 9173 


0. 2791 


0. 1297 


0. 0622 


4S 


2338. 03 


97. 4 


0. 8004 


0. 6991 


0. 4709 


0. 2577 


0. 1 187 


0. 0601 


49 


2397. 09 


99. 9 


0. 7210 


0. 6277 


0. 4289 


0. 2377 


0. Ills 


0. 0579 


SO 


2496. 07 


102. 3 


0. 6496 


0. 9667 


0. 3897 


0. 2190 


0. 1091 


0. 0996 



m 



STATIOM 4 


















K= 


62292. 


TEMP= 20 


. 00 a- 


493. 10 












AVG. 


VALUES: 1=1680. 13 


2- 17. 47 


Us Q. 017 


EX^ 141 


088 lETA- 


0. 0010 


RD=0. OOOC 


NO. 


TIMKHRS. 


,: DAYS) 


0.00 


0. 20 




0, 40 


0. 60 


0. 10 


I. 00 


1 


333. 60 


14. 9 


1. 1826 


0. 4855 


0. 


0350 


0. 0005 


0. 0000 


0. 0000 


1 


3S1. 63 


IS. 9 


1. 6873 


0. 7367 


0. 


0634 


0. 0012 


0. 0000 


0. 0000 


3 


407. 63 


17. 


2. 2709 


1. 0472 


0. 


1036 


0. 0025 


0. 0000 


0. 0000 


4 


433. 68 


18. 1 


2. 9130 


1. 4100 


0, 


1636 


0. 0048 


0. 0000 


0. 0000 


5 


459. 70 


19. 2 


3. 3903 


1. 9151 


0. 


2339 


0. 0016 


0. 0001 


0. 0000 


6 


483. 72 


20. 2 


4. 2809 


2. 2308 


0. 


3317 


0. 0142 


0. 0002 


0. 0000 


7 


311. 73 


21. 3 


4. 9635 


2. 7022 


0. 


4413 


0. 0223 


0. 0004 


0. 0000 


8 


337. 77 


22. 4 


5. 6210 


3. 1587 


0. 


3661 


0. 0332 


0. 0007 


0. 0000 


f 


363. 79 


23. 3 


6. 2397 


3. 6016 


0. 


7036 


0. 0473 


0. 0011 


0. 0000 


10 


389. 82 


S4. 6 


6. 8092 


4. 0422 


0. 


8513 


0. 0648 


0. 0018 


0. 0000 


11 


613. 84 


25. 7 


7. 3227 


4. 4319 




0060 


0. 0938 


0. 0021 


0. 0001 


12 


641. 86 


26. 7 


7. 7762 


4. 8317 




1646 


0. 1103 


0. 0042 


0. 0001 


13 


667. 89 


27. 9 


9. 1678 


3. 1775 




3243 


0. 1382 


0. 0060 


0. 0002 


14 


693. 91 


28. 9 


8. 4979 


3. 4867 




4823 


0. 1691 


0. 0083 


0. 0003 


IS 


719. 94 


30. 


9. 7679 


5. 7379 




6362 


0. 2029 


0. 01 11 


0, 0004 


16 


743. 96 


31. 1 


8. 9803 


3. 9906 




7840 


0. 2389 


0. 0145 


0. 0006 


17 


771.98 


32. 2 


9. 1390 


6. 1855 




9240 


0. 2768 


0. 0186 


0. 0009 


18 


791. 01 


33. 3 


9. 2472 


6. 3431 


2. 


0347 


0, 3160 


0. 0233 


0. 0012 


19 


824. 03 


34. 3 


9. 3092 


6. 4671 


2. 1732 


0. 3362 


0. 0216 


0. 0017 


SO 


882. 32 


36. 8 


9.3015 


6. 6274 


2. 


4054 


0. 4469 


0. 0427 


0. 0032 


21 


941. 02 


39. 2 


9. 1271 


6. 6475 


2. 


3772 


0. 3345 


0. 0596 


0. 0054 


m2 


999. 31 


41. 6 


B. 8302 


6. 3345 


2. 


6919 


0. 6148 


0, 0784 


0. 0084 


23 


lose. 01 


44. 1 


S. 4443 


6. 3742 


2. 


7541 


0. 6849 


0. 0914 


0. 0122 


24 


1116. 50 


46. 3 


S. 0004 


6. 1294 


2. 


7701 


0. 7430 


0. 1 188 


0. 0169 


2i 


1173. 00 


49. 


7. 5217 


3. 9393 


2. 


7472 


0. 7883 


0. 1387 


0. 0222 


S6 


1233. 49 


51. 4 


7. 0266 


3. 3200 


2. 6922 


0. 8210 


0. 1575 


0: 02Si 


27 


1291. 99 


33. 9 


6. 3292 


3. 1844 


2. 


6120 


0. 8415 


0. 1746 


0. 0344 


28 


1330. 48 


56. 3 


6. 0397 


4. 8423 


2. 


3126 


0. 8508 


0. 1896 


0. 0409 


2f 


1408.97 


38. 7 


5. 5657 


4. 3017 


a. 


3993 


0. 8502 


0. 2023 


0. 0473 


WQ 


1467. 47 


61. I 


3. 1122 


4. 1683 


2. 


2767 


0. 8410 


0. 2125 


0. 0533 


m 


1323. 96 


63. 6 


4. 6826 


3. 9462 


2. 


1495 


0. 8246 


0. 2203 


0. 0394 


32 


1584. 46 


66. 


4. 2787 


3. 5313 


2. 


0178 


0. 8022 


0. 2256 


0. 0648 


33 


1642:. 93 


68. 5 


3. 9015 


3. 2466 




8171 


0. 7750 


0. 2287 


0. 0696 


34 


1701. 43 


70. 9 


3. 5310 


2. 9721 




7514 


0. 7443 


0. 2296 


0. 0737 


33 


1739. 94 


73. 3 ■ 


3. 2269 


2. 7134 




6331 


0. 71 10 


0. 2286 


0. 0771 


36 


1818. 44 


75. 8 


2. 9283 


2. 4764 




5123 


0. 6751 


0. 2260 


0. 0798 


37 


1976. 93 


78. 2 


2. 6341 


2. 2350 




3969 


0. 6397 


0. 2219 


0. 0817 


38 


1933. 43 


SO. 6 


2. 4029 


2. 0305 




2873 


0, 6032 


0. 2165 


0. 0130 


3f 


1993. 92 


93. 1 


2. 1734 


1 . 9622 




1131 


0. 3668 


0. 2101 


0. 0836 


40 


2052. 41 


83. 5 


1 . 9642 


1 . 6894 




0866 


0. 3909 


0. 2030 


0. 0836 


41 


2110. 91 


88. 


1. 7737 


I. 3311 





9938 


0. 4939 


0. 1931 


0. 0831 


42 


2169. 40 


90. 4 


1. 6007 


I. 3864 





9111 


0. 4620 


0. 1869 


0. 0120 


43 


2227. 90 


92. 9 


1 . 4436 


1 . 2343 





8325 


0. 4294 


0. 1783 


0. 0105 


44 


2286. 39 


95. 3 


1. 3013 


1. 1341 





7598 


0. 3982 


0. 1695 


0. 0786 


4:9 


2344. 89 


97. 7 


1. 1724 


1. 0247 





6926 


0, 3686 


0. 1606 


0. 0765 


46 


2403. 38 


100. 1 


1 . 0331 


0. 9253 





6307 


0. 3406 


0. ISIS 


0. 0740 


47 


2461. 88 


102. 6 


0. 9S04 


0. 9331 





3738 


0. 3142 


0, 1430 


0. 0713 


48 


2520. 37 


103. 


0. 8333 


0. 7533 





3216 


0. 2894 


0. 1343 


0. 0613 


m 


2578. 87 


107, 3 


0. 7694 


0. 6793 





4737 


0. 2662 


0. 1261 


0. 0636 


50 


2637. 36 


109. 9 


0. 6920 


0. 6123 





4300 


0. 2446 


0. 1180 


0. 0626 



ii 
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